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Abstract 

A three-parameter discrete distribution is developed to describe 
the multiplicity distributions observed in total- and limited phase 
space volumes in different collision processes. The probability law 
is obtained by the Poisson transform of the KNO scaling function de- 
rived in Polyakov's similarity hypothesis for strong interactions as well 
as in perturbative QCD, ^(z) oc z"exp(— z^). Various characteristics 
of the newly proposed distribution are investigated e.g. its generating 
function, factorial moments, factorial cumulants. Several limiting and 
special cases are discussed. A comparison is made to the multiplicity 
data available in e"^e~ annihilations at the peak. 
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1. Introduction 



Over the past years much attention has been focused on the description of 
multiphcity distributions (MDs) for hard processes in perturbative quantum 
chromodynamics (pQCD), see ref. [1] for recent reviews. Already in the 
lowest double logarithmic approximation (DLA) it was demonstrated that 
QCD provides a natural explanation of the KNO scaling law 

Pn=^^^{^^ (1-1) 

(n) \{n)) 

of the MDs [2] but the width of the scahng function t/j(z — n/n) was overesti- 
mated. To improve on the agreement with experiments higher-order pertur- 
bative corrections should be taken into account. Mathematically this mani- 
fests itself in the appearance of a new expansion parameter 7g, the product 
of the QCD multiplicity anomalous dimension 7 oc y/aj and the rank q of the 
moments of P„. For example, the influence of recoil effects [3] on multiplicity 
fluctuations is seen most clearly when one considers the normalized factorial 
moments 

(n(n - l)...(n-g + 1)) 
■ (1-2) 

In the next-to-next-to-leading log approximation (NNLA) the Fq are shown [4] 
to be dominated for large ranks q by the F-function of the rescaled rank qf/i, 

Fq^T{\ + q/ji)D-'i (1.3) 

where = (1 — 7)"^ and D is a scale parameter depending on /x. Since 
7 ~ 0.4 at the peak, the strength of higher-order multiplicity fluctuations 
are drastically reduced by recoil effects and therefore the tail of the KNO 
function gets modified. Instead of the exponential fall-off predicted in the 
DLA, a faster than exponential decay law emerges for large z [4] 

i^iz) (X z"" ex^{-[DzY) (1.4) 

with a — ?>ijl/2 — 1 and ~ 5/3. This type of behaviour reduces the width 
of ip{z) overestimated by the early DLA calculations. 

Another manifestation of higher-order pQCD effects on the MDs is best 
seen in the g-dependence of the ratio Hq = Kg/Fg invented by Dremin [5] 
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where Kq denotes the normahzed factorial cumulant moments, 

Kq^Fq-Y,^i^~^^Kq_,F,. (1.5) 

In the DLA Hq decreases monotonously as q~'^ but the inclusion of higher 
order perturbative corrections yields quite nontrivial ^'-dependence. In the 
next-to-leading log approximation (MLLA) Hq acquires a negative minimum 
at g ~ 5 and approaches for large ranks q the abscissa from below. In 
the NNLA the negative minimum at g ~ 5 is followed by sign-changing 
oscillations of the moment ratios [6] . Quite recently this peculiar behaviour 
has been experimentally confirmed in e^e~ annihilations [7]. 

In a historical context it is worth mentioning that some of the above pre- 
dictions of pQCD were put forward already in 1970 by Polyakov formulating 
a similarity hypothesis for strong interactions in e^e~ annihilations [8]. He 
obtained asymptotic KNO scaling behaviour for the MDs (two years earher 
than the original KNO-paper) with a scaling function that behaves according 
to 

rn^lT^'"'""'''^ ""^n" (1-6) 
^ ' as z — s> ^ ' 

where < 5 < 1/2 and the small- 2; behaviour of the scaling function was 
found, with certain assumptions, to be a monomial [8]. Clearly, Polyakov 
arrived in his model at the same law for multiplicity fluctuations as the pQCD 
result (1.4). A KNO scaling function of the form of Eq. (1.4) was studied in 
more or less detail by others as well [9] most notably by Krasznovszky and 
Wagner. We call attention also to the remarks made by Koba, Nielsen and 
Olesen in connection with the moment problem [2] to which we shall return 
later. 

The goal of the present paper is the development of a discrete distribution 
well suited to describe available multiplicity data when the popular negative 
binomial parametrization fails. Our guiding principle is the incorporation of 
the pQCD-based characteristics i) factorial moments Fq dominated by the 
F-function of the rescaled rank q/ fx, ii) factorial cumulant-to- moment ratios 
Hq displaying nontrivial sign-changing oscillations, and Hi) distributions P„ 
possessing squeezed high- multiplicity tail in accordance with Eq. (1.4). We 
shall see that a simple generalization of the negative binomial law with one 
additional parameter satisfies all these requirements. 
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2. The scaling function in pQCD 



Let us first consider in more detail the modified KNO scaling function (1.4) 
obtained in pQCD with fixed coupling approximation [4]. For later conve- 
nience we use a somewhat different choice of the parameters. Introducing 
the shape parameter k = {a + l)//! and denoting the scale parameter by A 
we get 

ij{z) =Mzi'^-^ exp {-[\zY) (2.1) 

where A/" is a normalization constant. The scaling function tlj{z) should satisfy 
the two normalization conditions 

/ ij(z)dz= zij(z)dz=l. (2.2) 
Jo Jo 

The first condition determines the constant jV. It can be obtained by the 
use of the integral [10] 



fz'^''-^e.p{-[Xzr)dz^^ (2.3) 



providing the reciprocal of A/". Thus the complete form of the scaling function 
(at the moment without the second normalization constraint) is 

m = ^ >^''^""' exp (-[Xz^) . (2.4) 

Eq. (2.4) is known in the mathematical literature as the generalized gamma 
distribution [11]. Obviously, the fi = 1 special case yields the ordinary 
gamma distribution which corresponds in pQCD to the KNO scaling function 
at asymptotically high energies (7 oc ^/a^ — > 0) . 

The probability density function (2.4) has a long history [12]. Its first 
appearance in 1925 is due to L. Amoroso who analyzed the distribution of 
economic income. The next application is related to the grinding of materials: 
in 1933 Rosin, Rammler and Sperling arrived at Eq. (2.4) as the size distri- 
bution of grains produced by comminution. Nevertheless in the literature of 
particle size distributions the formula is named after Nukiyama and Tana- 
sawa who rediscovered it in 1939 studying drop size distributions in sprays. 
Widespread use of Eq. (2.4) by mathematicians started only in 1962 by the 
paper of Stacy. He introduced it as the generalized gamma distribution and 
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since then this name is the most frequently used one. Some of the recent 
apphcations include polymer physics [13] and the theory of multiplicative 
processes [14]. 

The moments of the scaling function (2.4) are readily obtained from the 
integral given by Eq. (2.3) which yields 

{z^) = J^ .V(.)d.= ^^-. (2.5) 

Eq. (2.5) with q = 1 and the second normalization condition in Eq. (2.2) 
restrict the scale parameter to 

which completes the analytic form of the scaling function. Let us now 
consider the cumulative distribution corresponding to the probability den- 
sity (2.4). It is given by 

where 7(-) denotes the incomplete gamma function. Eq. (2.7) is particularly 
useful in connection with KNO-G scaling [15,16] i.e. if the discrete mul- 
tiplicity distributions P„ are approximated by the continuous scaling func- 
tion ip{z), 

iIj{z) dz = ^{Zn+l) - (p{Zn)- (2.8) 



Although KNO-G scaling is becoming increasingly popular in the description 
of MDs, with ip{z) chosen to be a shifted lognormal density, we will follow a 
different approach to construct P„ from the scaling function. 

Perhaps the most advantageous property of the generalized gamma dis- 
tribution is that ■i/j{z; k, A, /i) can be used to specify a variety of well known 
probability laws. Some examples are cited in Table 1 with the corresponding 
set of parameters. It is seen that the exponent /i can take negative values as 
well in which case the normalization constant in Eq. (2.4) involves [//[. More- 
over, for negative /i the moments {z'^) of the distribution are finite only for 
ranks q satisfying q/ /J. > —k. A good example is the completely asymmetric 
Levy law of index a = 1/2 having no finite, positive-rank moments. Finally 
let us call special attention to the last two rows of Table 1 that display the 
limit distributions of ip{z; k, A, /i) involving for — > the lognormal law. 
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3. MDs at preasymptotic energies 



Comparing the KNO function (2.4) to the pQCD result quoted in Section 1 
we see that the more accurate account of recoil effects in [4] yields for ip{z) 
a generalized gamma distribution with fixed shape parameter k — 3/2. Fur- 
thermore, the factorial moments given by Eq. (1.3) are the same as the 
ordinary moments (2.5) of ip{z; |, A, /x) which is due to the identification of 
{z'^) and Fq in ref. [4]. Their equivalence holds valid for the /i — 1 asymptotic 
energy scale as well as for preasymptotic energies characterized by // > 1 if 
Pn is defined by the Poisson transform of tlj{z), 

ij{z)^—^e-'''dz. (3.1) 

77/. 

For Eq. (3.1) the asymptotic scaling limit of the MDs is [17] 

hm P„ - y\ ^ (^] . (3.2) 

n/(n) fixed ^ ' \ ^ ' / 

Thus the knowledge of the shape of P„ at preasymptotic energies enables one 
to guess the asymptotic shape of ip{z). We mention that the scaling limit (3.2) 
naturally arises in Polyakov's model [18] suggesting the use of (3.1) instead 
of (2.8) to construct P„ from ^p{z). 

The probabihty generating function of P„ defined by the Poisson trans- 
form (3.1) reads as follows: 

g{u) = Y.{l-uYPn= iPiz)e-''''''dz, (3.3) 

i.e. g{u) is given by the Laplace transform of ip{z) [17]. Unfortunately the 
Laplace transform of Eq. (2.4) can be expressed in terms of special functions 
only for a few specific values of the parameter /i. To go further we have to 
make use of Fox's generalized hypergeometric function, H(a;). The reader 
unfamiliar with the theory of Fox functions can find a compilation of the 
necessary formulae in the Appendix. The KNO scaling function (2.4) is 
expressed in terms of H(x) as the following H-function distribution: 



m - Hi? 



Xz 



{k - 1//X, 1//.) 



(3.4) 
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The above form is obtained from Eq. (A. 7) using the identity (A. 4) and the 
integral (2.3) to ensure proper normahzation, see also [22]. The probability 
generating function of P„ defined by the Poisson transform of Eq. (3.4) can 
be easily evaluated with the help of (A. 5), the Laplace transform of H{x). It 
yields 



and 



g{u) 



g{u) 



m 



A 



un 



un 

T 



(1, 1) 



< < 1 



(1 - k, llix) 
(0, 1) 



\x>\ 



(3.5) 



(3.6) 



where A is given by Eq. (2.6). The necessity of two separate expressions for 
Q{u) follows from the existence conditions of H(a;) discussed in the Appendix. 
The /i = 1 case is of course the generating function of the negative binomial 
distribution, Q{u) = {l+un/k)~^. The < case will be studied elsewhere. 

The characteristic function (f){t) of the generahzed gamma density (2.4) is 
usually written as an infinite sum [21]. In terms of the Fox function (f){t) can 
be expressed in a much simpler way by changing the variable in Q [u) given 
by Eqs. (3.5-6), 



POO 

/ ^izy^'dz = g 

Jo 



un 



-it) 



(3.7) 



with, in general, unconstrained scale parameter A. For n > 1 (f){t) was 
obtained earlier in ref. [22]. The probability generating function provides 
also the Poisson transform of tp{z) for n = through Pq = 

Let us now consider the Poisson transform of (3.4) for arbitrary n. By 
simple algebra, using identities (A. 2-4) and the Laplace transform (A. 5) of 
H(x) one arrives at the following discrete H-f unction distributions: 



P 

" n\r{k) 



H};1 



(l-n, 1) 
(k, 1//X) 



, < // < 1 



(3.8) 



and 



" n\r{k) 



;i - k, i/f,) 

{n, 1) 



, li>l. 



(3.9) 



The A* = 1 case yields the Poisson transform of the ordinary gamma distri- 
bution, i.e. the negative binomial law. The splitted parameter space for /i 
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reflects an important difference between the two expressions for P„. In case 
of Eq. (3.8) Pn is infinitely divisible, just as the negative binomial for = 1, 
whereas the // > 1 case given by Eq. (3.9) violates this feature. Due to 
the preservation of infinite divisibility under Poisson transforms the same 
distinction between the two //-domains holds for the generalized gamma dis- 
tribution [23]. 

The factorial moments of P„ can be determined through the equivalence 
of {z'^) and Fq for Poisson transforms [17]. Plugging into Eq. (2.5) with scale 
parameter A restricted by (2.6) we get 

r(^+?/^ 1 

As intended, the factorial moments are dominated by the F-function of the 
rescaled rank for large q. With the help of Eq. (1.5) we have calculated 
the moment ratios Hq = Kq/Fq for shape parameter k = 3/2. The behaviour 
of log \ Hq\ over the ^-q plane is shown in Fig. la. The peculiar g-dependence 
for // > 1 is due to sign-changing oscillations of Hg. Observe that the pattern 
of oscillations is nontrivial, i.e. not alternating as q takes even/odd values, 
see Fig. lb. Qualitatively similar sign-changing oscillations of Hq occur for 
a different choice of the shape parameter k. 



4. Comparison to e+e data 

One of the advantages of P„ given by Eqs. (3.8-9) is its generality: the Poisson 
transform of many well known probability laws (e.g. of those cited in Table 1) 
is a special case. The negative binomial distribution (/x = 1) not rarely fails to 
give a reasonable description of the MDs. For example, in e+e~ annihilations 
the full phase-space MDs at the peak show significant deviation from 
a negative binomial shape as reflected by the x^/d.o.f. = 66/23 and 68/24 
values reported by the Delphi and SLD collaborations [24,7]. Moreover, the 
sign-changing oscillations displayed by the SLD data for Hg [7] can not be 
reproduced by a (possibly truncated) negative binomial. It is natural to ask 
how the fjLy^l special cases of Eqs. (3.8-9) can describe the observed features. 

To answer this question we carried out fits to the e~^e~ multiplicity data 
at the Z'^ peak [7,24]. In the fitting procedure numerical evaluation of the in- 
tegral (3.1) is implemented using 96-point gaussian quadrature. The scaling 
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function (2.4) is log-linear so that y = lnz can be written as y = w/fi — \nX 
where w has probability density f{w) = exp{kw — e'^)/T{k). Making the pa- 
rameter transformations p — k~^^'^, a — p/ ji and a = In A, further, reflecting 
the model about the origin p = to negative p the scaling function (2.4) 
takes the form 

\v\ r 1 

Mz) = exp p-^w -e"" if p 7^ (4.1) 

with w — p{hiz -\- a)/(T. The reparametrization allows the lognormal law to 
be mapped to the origin, 



, , , . 1 (ln^ + 

V'(^) = exp ■ 

'27: (TZ 2: a 



2 



if p - 0. (4.2) 



Above, the location parameter a is restricted by the second normalization 
condition in Eq. (2.2) to 

^ ^ r In r(p-2 + a/p) - In r(p-2) if p ^ 

\ 0-2/2 if p = ■ ^ ' ' 

The distributions represented by Eqs. (4.1-2) are stochastically continuous 
in the parameters and include besides the lognormal (p = 0) e.g. the gamma 
(p = cr), exponential (p = cr = 1) and WeibuU (p = 1) densities. For further 
details on the above transformation of generalized gamma variates the reader 
is referred to [25] . 

In the numerical evaluation of the integral (3.1) we have used Eqs. (4.1-2) 
for 'i\}{^z). There is a strong reason of replacing the parametrization (2.4) 
with those of (4.1-2) in the fitting procedure. The original fits using the 
form (2.4) were able to reduce significantly the corresponding to the = 1 
negative binomial case but they resulted unexpectedly small values of For 
example, fitting the Delphi data yields 'x^ j^^.o.i. = 33/23, k ^ 1000 and 
/i ^ 0.1 indicating that the observed P„ is a Poisson transformed lognormal 
distribution. Performing the fits using (4.1-2) reinforced our finding for each 
data set. This can be seen from Table 2 where the outcome of the fits are 
collected. In Fig. 2a the best-fit theoretical P„ is displayed for the SLD data. 

The above result seems to be in sharp contradiction with the SLD data 
for the moment ratios Hq [7]. Sign-changing oscillations can occur in our 
model only for > 1, if P„ is not infinitely divisible. Fitting the Hq data 
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with the help of Eqs. (3.10) and (1.5) produced unacceptable x^. To allow 
the possible influence of truncation effects [26] we recalculated the factorial 
moments Fq in terms of P„ using = 6 and Umax — 50 in accordance 
with the SLD data [24]. The resulting behaviour of Hq is shown in Fig. 2b, 
the agreement between the theoretical curve and the data points is excellent. 
The corresponding and the values of the best-fit parameters arc cited in 
the last row of Table 2. Clearly, the obtained parameters are very similar 
to those found in the previous fits. Although there is a small but significant 
deviation fromp = we can safely conclude that the multiplicity data in e+e~ 
annihilations at the peak favour a lognormally shaped KNO function 'il}{z) 
in the asymptotic limit given by Eq. (3.2). 

5. Summary and conclusions 

We have developed a three-parameter discrete distribution to analyse the ex- 
perimental data for MDs in different collision processes. It is obtained by the 
Poisson transform of the generalized gamma density (2.4) and thus it provides 
a generalization of the popular negative binomial distribution. The analytic 
form of Pn is expressed in terms of H-functions. It is given by Eq. (3.8) for 
infinitely divisible /-"„ and by Eq. (3.9) if P„ does not obey this feature. The 
11=1 marginal case between the two families of distributions is the negative 
binomial law. Violation of infinite divisibility for Eq. (3.9) allows nontrivial 
sign-changing oscillations of the factorial cumulant-to-moment ratios Hq. For 
> 1 we thus have the ability to reproduce this peculiar g-dcpcndcncc of Hq 
data without truncation. Since the factorial moments (3.10) arc dominated 
by the F- function of the rescaled rank q/ ^ for large q we can also repro- 
duce possible enhancement (// < 1) and suppression (/x > 1) of multiplicity 
fluctuations with respect to the negative binomial behaviour. 

Fitting the newly developed distribution to the e"'"e~ multiplicity data 
at the peak we have found /i < 1 departure from the negative binomial 
law. In fact the observed MDs are in agreement with the Poisson transform 
of the /X ^ limit distribution of (2.4) which is the lognormal density. 
Even the sign-changing oscillations of the SLD data for Hq can be fully 
reproduced after taking into account truncation effects. It is important to 
emphasize that the ^ limit manifests itself through the factorial moments 
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of Pn being equivalent to the ordinary moments of the asymptotic ip{z). At 
the peak the KNO function exhibits suppressed high-multiphcity tail in 
accordance with the pQCD results quoted earlier. We disagree with ref. [16] 
stating that ip{z) is lognormally shaped at current energies, which became 
a common belief in recent years. Fitting the scaling function (2.4) to ip{z) 
observed by experiments displays clear deviation from lognormality yielding 
/c ~ 5 and ~ 1.6 best-fit parameters. The marked difference between the 
experimental and asymptotic scaling functions warns that the ordinary and 
factorial moments of P„ can not be identified at current energies and sign- 
changing oscillations caused by pQCD effects are more likely to occur for the 
ordinary cumulant-to-momcnt ratios. 

The Poisson transformed lognormal shape of MDs calls attention to the 
moment problem considered already in ref. [2] for Eq. (2.4). According to 
our results the KNO function at the asymptotics (3.2) is not determined by 
its moments, a well known feature of the lognormal law. There is another 
important aspect of lognormality at asymptotic energies. It suggests that the 
limiting shape of ijj{z) can be guessed on more elementary grounds than the 
perturbative treatment of parton evolution equations. The multiplicative 
nature of parton cascades in QCD and the large number of steps of the 
cascade processes at very high energies cause the central limit theorem to 
come into operation determining uniquely the asymptotic shape of the scaling 
function. 
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Appendix 



Here we give a brief summary of the basic properties of Fox's generalized 
hypergeometric function, H(a;). It heavily relies on article [19] and on the 
book [20] where the interested reader can find more details. 

The H-function of Fox is defined in terms of a Mellin-Barnes type integral 
as follows: 



X 



\Clpi Oipj 

(b,, 13,) 



^/^AW.-d. (A.1) 



where x ^ and 

x^ — exp {s [In \x\ + i arg(x)]} 
in which arg(x) is not necessarily the principal value. Further, 

nf=i rjbj - f3js) UU r(l - aj + ajs) 



nUm+i r(i - bj + M nU+i r(% - ajs) 

where m, n,p, q are integers satisfying 

0<n<j9, 1 < m < q. 



The parameters (ai, . . . , a^) and . . . 
and (/?!,..., Pq) are positive numbers, 
unity. The parameters are restricted 



bq) arc complex, whereas {ai, . . . , ctp) 



An empty product is interpreted as 
by the condition that aj{bh + v) 
Phidj — 1 — A) for I/, A = 0, 1, . . .; h — 1, . . . , m; j — 1, . . . ,n. The con- 
that the points s = {bh + i^) / Ph and 
left of C respectively while C extends 
A; is a constant with k > |Im bh\/ Ph- 
the following two existence conditions 



tour C in the complex s plane is such 
s = {ttj — 1 — u) / aj lie to the right and 
from s = oo — ik to s = oo + ik where 
The H-function makes sense only if 
are satisfied: 

i) X and p > where 



P 



ii) p — and < |a;| < ^ ^ where 



j=i i=i 
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Under these conditions H(a;) is an analytic function for a; 7^ 0, in general 
multivalued, one- valued on the Riemann surface of Inx. 

Elementary properties of H(x) utilized in the body of the paper: 



UJm,n 



c |_|m,n 



p, q 



X 



X 



x"- 





P,) _ 


UJn,m 
= "qlp 


1 

X 


(1 - b,, Pg) 
(1 - ap, ap) _ 




ap) 


p> q 


X 


(cip -|- C Op, Op^ 


K 


P,) _ 


{bq + C(3q, (3q) 


{dp, 

{K 


ap) 




q 


^ (cp, ap/c) 

{bg, Pg/c) _ 



(A.2) 

(A.3) 
(A.4) 



where c > 0. The first identity, Eq. (A.2), is an important one because it 
enables us to transform H(a;) with p < to H(l/x) with p > and satisfying 
the existence condition i) given above. The Laplace transform of H(a;) reads 
as follows: 



r 

Jo 



^P,q 



( cx) e 



1 

c 



-rxj _ Jll |_|n+l,m 
— "q, p+1 



1 



(0, 1) 



- bg - Pg, Pg) 

'1 — a 



ap , ap j 



(A.5) 



Special cases of the H-function include such as Meijer's G-function, Bessel, 
Legendre, Whittaker, Struve functions, the generalized hypergeometric func- 
tion and several others. Meijer's G-function is obtained when the H-function 
parameters (ai, . . . , ap) and . . . , f3q) are unity. The generalized hyperge- 
ometric function pFq is related to H(a;) through 



pFq{ctp, bq, x) — 



ul,P 



-X 



(1 - Op, 1) 
(0, 1), {1-hq, 1) 



(A.6) 



Prom this relationship one can obtain numerous distributions of statistics 
(associated e.g. with the Gauss- and confluent hypergeometric functions 2-^1 
and iFi) in terms of V\{x). Two further special cases of H(x) frequently 
encountered in statistics are 



and 



x°exp(-Ax^) = A-"/^ hJ'? 



x"/(l + Xx^) = A-"/'' 



Xx^ 


HP, 1) . 


Xx^ 


{a/(3, 1) " 
HP, 1) . 



(A.7) 



(A.8) 
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For a rich collection of particular cases of the H-function, see ref. [20]. 

Let us now consider the continuous random variable x G (0, oo) whose 
probability density function is given by 



\x 



(A.9) 



with scale parameter A > and normalization condition /q°° f{x)dx = l. The 
random variable x is called a H-function variate and the probability law (A.9) 
with f{x) > is the so-called H-function distribution. Its parameters (op, ap) 
and {bq, /3g) should satisfy all the restrictions given earlier in this Appendix. 
The characteristic function of f{x) can be obtained from Eq. (3.7) and the 
Laplace transform (A. 5): 



it 



(0, 1), (1 - Op - ap, ap) 



The gth moment of the H-function distribution is related to its r 
Mellin transform through 



ra-i 



M nr=ir(&,+r/j,)n-=ir(i-a,- 

nl=n.+i r(l - - r(5^) Y{%n^i r(% + ra^) ' 



(A.IO) 



(g+l)th 



(A.ll) 



Many of the classical probabihty densities are special cases of the H-function 
distribution and can be written in the form of (A.9). This was utilized in 
Eq. (3.4) which enabled us to express the Poisson transform of the scaling 
function (2.4) analitically as a discrete H-function distribution. 
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FIGURE CAPTIONS 

Fig. la: Sign-changing oscillations of Hq = Kq/Fq with Fq given by Eq. (3.10) 
for k — 3/2. The neighbouring bumps at a fixed q are ifq-intervals of oppo- 
site sign in yU. For ^ < 1 Hq is always positive. The yU-scale is logarithmic 
and for clarity only the odd-rank moment ratios are shown. 

Fig. lb: Shces through Fig. la with q = 5, 20 (left) and with = 1, 3/2 
(right). 

Fig. 2a-b: The best fit to the SLD data for P„ and Hq as described in the 
text. The fit parameters and the quality of fits are shown in the last two 
rows of Table 2. 
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Distribution 


shape par. 


scale par. 


exponent 


Generalized gamma 


k 


A 




Gamma 


k 


A 


1 


Chi-square, n d.o.f. 


n/2 


1/2 




Exponential 


1 


A 




WeibuU 


1 


A 


1^ 


Stratonovich 


k 


A 


2 


Rayleigli 


1 






Maxwell molecular speed 


3/2 






Maxwell mol. velocity 


1/2 






Chi, n d.o.f. 


n/2 


1/V2 




Half-normal 


1/2 






Circular normal 


1 






Spherical normal 


3/2 






Pearson type V 


k 


A 


-1 


Asymmetric Levy, a = 1/2 


1/2 






Lognormal 


k ^ oo 


A 




Pareto 


k 


A 





Table 1. Special- and limiting cases of the generalized gamma density (2.4) 



EXPEIIIMENT 


vVd.o.f. 


l> 


(") 


a 


Aleph 


4.2/24 


-0.230 ±0.311 


20.970 ±0.407 


0.200 ±0.014 


Delphi 


30.2/23 


-0.089 ±0.113 


21.311 ±0.140 


0.199 ±0.006 


L3 


14.9/20 


-0.123 ±0.238 


20.660 ±0.349 


0.205 ±0.016 


Opal 


10.7/24 


-0.124 ±0.130 


21.370 ±0.128 


0.208 ± 0.008 


SLD 


31.2/20 


-0.001 ± 0.060 


20.892 ± 0.068 


0.204 ± 0.003 


SLD (Hg) 


16.6/13 


0.060 ± 0.029 


21.346 ±0.126 


0.214 ± 0.004 



Table 2. The results of fits to the e+e multiplicity data at the peak 




Figure la 





Figure 2b 



